Stripe Formation within SO(5) Theory 
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We study the formation of stripe order within the SO (5) theory of high-T c superconductivity. Spin 
and charge modulations arise as a result of the competition between a local tendency to phase 
separate and the long-range Coulomb interaction. This frustrated phase separation leads to hole- 
rich and hole-poor regions which are respectively superconducting and antiferromagnetic. A rich 
variety of microstructures ranging from droplet and striped to inverted-droplet phases are stabilized, 
depending on the charge carrier concentration. We show that the SO(5) energy functional favors 
non-topological stripes. 
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One of the most striking features of the cuprates is the 
proximity between antiferromagnetic (AF) and supercon- 
ducting (SC) phases as a function of doping. Recently it 
has been proposed that these two phases are unified by 
an approximate SO(5) symmetry jlj. A number of ex- 
perimental consequences of this theory have been worked 
out (§]-[§■ Although SO (5) appears to be a natural frame- 
work for understanding the cuprates, no experiment has 
unequivocally tested the fundamental validity of the the- 
ory. One of its most direct predictions is the existence 
of a first-order transition from the AF to SC state as the 
chemical potential fi is increased beyond a critical value. 
However, this prediction is complicated by the fact that 
the doping x (not p) is the experimentally tunable pa- 
rameter. Experimentally, it is found that in the vicinity 
of the AF/SC transition, the cuprates show an increased 
sensitivity to disorder and inhomogeneity. In this Letter, 
we study this region of the phase diagram in the presence 
of the long-range Coulomb interaction within the SO(5) 
formalism and show how spatially inhomogeneous states 
can emerge. 

In the T — (i phase diagram of SO (5) theory, there is a 
first-order line separating the AF and SC phases, across 
which the charge carrier density x jumps discontinuously. 
In the T — x phase diagram, this translates into a two- 
phase region where AF and SC phases coexist. Phase 
separation into hole-rich and hole-poor regions was also 
noticed in studies of the t — J model However, as 

Emery and Kivelson M argued rather successfully, the 
long-range Coulomb interaction between charge carriers 
prevents macroscopic phase separation. The competi- 
tion between the local tendency toward phase separation 
and the long-range Coulomb interaction leads to modu- 
lated domain structures at mesoscopic scales ]To|— |l3[| . In 
the SO (5) theory, the hole-rich and hole-poor regions are 
respectively identified as having a superconducting and 
antiferromagnetic character. The spin and charge mod- 
ulations of the system are interpreted as textures of the 
SO (5) superspin as it rotates in SO (5) space. 

There is considerable evidence for modulated mi- 



crostructure in the oxides. Domain formation has 
been reported in La2- a; Sr a ;Cu04 (LSCO) in muon 
spin resonance jL4j, NMR and neutron diffraction ex- 
periments Jl~5| ]. Neutron scattering measurements in 
Lai.6-xNdo.4Sra;Cu04 (LNSCO) provide direct evidence 
for stripe ordering in which the phase of the AF order 
shifts by 7r across a domain wall . Furthermore, recent 
inelastic neutron scattering measurements in underdoped 
YBasCuaO^ (YBCO) Q and ARPES measurements 
in underdoped Bi 2 Sr 2 CaCu20 8 _ ;E (BSCCO) @ are not 
inconsistent with a striped phase interpretation. 

In the mean field approximation of the SO (5) theory 
one minimizes the classical energy 



H\{n a ,p a ) = - 



ab 



E r 



(1) 



with L a b = n a p — n p a and constraints n a — l, n a p a — 
0. L a b and n refer respectively to the SO(5) gener- 
ators of rotation and the 5-component superspin n = 
(Re(A),N x ,N y ,N z ,Im(A)) B). The last term E c is the 
Coulomb energy. Since, in SO (5) theory, the hole density 
is given by L15, the charge density is p(r) = £15 (r) — ex, 
where ex is the charge of the neutralizing counterion 
charges which are assumed to be static and homoge- 
neously distributed with a density x. We conjecture then 
that, if SO(5) theory is still valid in the presence of the 
long-range Coulomb interaction, /citecoul the Coulomb 
energy in the mean field approximation will be given by 

E c = ^JJp(r)V c (r-r')p(r')drdr' . (2) 

It is important to emphasize that for homogeneous 
phases one recovers the basic SO(5) model because the 
charge density vanishes exactly making the Coulomb in- 
teraction irrelevant. Hence the influence of the Coulomb 
term arises solely in the phase separation regime. Sec- 
ond, the assumption of immobile static counterions is 
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known to fail for La2Cu04+a p0|| . In this case, the oxy- 
gen ions are mobile enough as to screen the charge inho- 
mogeneities which leads to macroscopic phase separation 
of superconducting and antiferromagnetic domains. This 
last fact provides strong evidence for the correctness of 
the SO (5) picture. 

It can be proved that if H\ is symmetric with respect to 
rotations within AF and SC subspaces (xi5 — Xs,X23 — 

X34 = X24 = Xa, X12 = Xl3 = XU = X23 = X24 = Xtt), the 

minimal configuration has a form ft = (n\, n 2 , 0, 0, 0) = 
(cos 9, sin#, 0,0, 0), p — (0, 0,0, 0,p5) and the constraints 
are automatically satisfied if we use the variables (6,P5). 

With the addition of the long-range Coulomb interac- 
tion to the Hamiltonian, the behavior of the system is 
no longer tractable analytically and one must resort to 
numerical analysis. We note that a classical spin Hamil- 
tonian 



H 2 = jJ2Si-S j -2Kj2(Sf) 1 



(3) 



can be transformed into the local part of H\ by a Hal- 
dane map (see for details). Here J > and K > 
are known functions of Xab,9,Ps- To match the whole 
expression (uh we add a term 



2 ^ 



x)(S! - x) 



(4) 



where e is the dielectric constant of the material. 

As emphasized earlier, since experiments are per- 
formed at constant carrier concentration, the Hamilto- 
nian E = H 2 + V is subject to the doping constraint 
(S z ) — x. Numerically it is easier to study E than Hi 
because the hole density is given explicitly by S z rather 
than implicitly by L15. The properties of E are stud- 
ied using Monte Carlo simulations. In order to find the 
lowest energy state of the system, we perform simulated 
annealing from high-temperature. We assume an N x N 
2-dimcnsional lattice where N can be up to 40 unit cells. 
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FIG. 1. Phase diagram as a function of strength of the 
Coulomb interaction 1/e and the hole density x for K = 0.4 
and J = 1 which gives x c = 0.5. 



In the absence of the long-range Coulomb interaction, 
one can easily show that the system phase separates for 
densities x less than x c = K/ (2 J — K). The addition of 
the Coulomb term leads to a rich variety of modulated 
structures, which are shown on the phase diagram of Fig.. 
For large dielectric constant, three phases are found to be 
stabilized: a droplet phase made of SC droplets embed- 
ded in an AF background, a striped phase of alternating 
SC and AF stripes and an inverted-droplet phase where 
the droplets are antiferromagnetic. Our numerical solu- 
tions (Fig j^) show that the superspin stays in the AF or 
SC directions inside the domains and changes only in the 
thin domain walls. The structure represents a collection 
of solitons rather than a small modulation of the direc- 
tion of n. The superconducting density switches between 
and x c in AF and SC domains which means that the 
superconducting area fraction A sc /A — x/x c leading to 
a linear relation between the superfluid density and the 
doping as seen in some experiments 22 . A simple physi- 
cal argument for the pattern shape can be given in terms 
of interface energy [^3| . As long as one of the phases ( AF 
or SC) is in the minority, the energy of the AF/SC inter- 
face predominates over the Coulomb energy and circular 
domains are preferred as they minimize the length of this 
interface. However, the situation is reversed for x ~ x c /2 
where the repulsive interaction leads to dipole formation 
which favors elongated domains such as stripes. 

The striped phase is reminiscent of the domain struc- 
ture observed in LNSCO, though the rows of charge are 
superconducting in our model. It is interesting to con- 
sider the superspin texture in the striped phases, namely, 
the relative phase shift of every other stripe. Numerically, 
it is found that the lowest energy states do not show any 
winding of the superspin in space. Hence, the phase of 
the AF order parameter does not shift by tt on cross- 
ing a SC stripe. The same results were obtained for a 
simulation on a spin ladder where all configurations are 
necessarily one-dimensional. 

This absence of topological phase shift is in striking 
contrast to experimental data |]l6| ; however, it can be 
proved analytically for the minimal periodic 1-D config- 
uration of (p]) using the theorem of Pryadko et al. J|J] . 
Theorem: For a functional 

£ = J(^;) 2 +E loc (v 2 (r),r)dr 

+ Jj p(v 2 (r),r)V(r~r r )p(v 2 (r'),r')drdr' (5) 

of the function of one argument v(r) the minimal config- 
uration under a constraint j p(v 2 (r),r) dr = 0, does not 
cross zero. 
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FIG. 2. Geometric phases as a function of doping: Hole 
density profile for the (A) Droplet phase [x = 0.12), (B) 
Striped phase (x = 0.24) and (C) Inverted-Droplet phase 
x = 0.35. The parameters chosen are J = 1,K = 0.4 and 
e = 6. 

Applying it to (|l|) we note that as a function of 9 both 
p and H\ can be expanded in even powers around the 
points 9 = 0, 7r and thus 9 m in(r) can not cross these 
levels. In addition, let us perform a variable change 
(9,P5) — > (6,q = ps/cosfl). Its Jacobian is sometimes 
infinite, but not on the minimal solution for which mini- 
mization of Hi with respect to p§ gives: 



— cost 



P 5 



V c {r - r')p{r')dr' (6) 



cos 2 9/xs + sin 2 9/xir 

After the variable change Hi(9,q) can also be expanded 
in even powers around the points 9 = 7r/2,37r/2, so the 
minimal solution does not cross these levels either. Al- 
together, 9 m i n (r) always stays in one of the four quad- 
rants of the circle. For low dielectric constant, i.e. weak 
screening, we find that phase separation is precluded alto- 
gether. The system exhibits a homogeneous mixed phase 
in which the superspin points neither purely in the AF 
or SC direction. This state is reminiscent of the putative 
supersolid phase in 4 He as both order parameters (AF 
and SC) are nonzero everywhere in the sample. 
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chemical potential \i 

FIG. 3. x versus fi curve for three different values of e, for 
paramteters J = 1 and K = 0.4 which yield x c = 0.5. The 
curves e = (homogeneous) and e = oo (phase separation) 
are obtained from the analytic solution of Eq. . The circles 
are the data of the numerical solutions for e = 20. 

After the minimum of E is found, the chemical poten- 
tial can be numerically calculated as fi = dE/dx. As 
shown in Fig. ^|, p{x) becomes non-monotonic and a re- 
gion ofdpi/dx < appears. Such a region is prohibited in 
thermodynamics, but in models with continuous charge 
density (as opposed to point charges) it is generic. The 
origin of the region in which the chemical potential is 
double- valued is illustrated by Fig.[|. In the absence of 
Coulomb interactions and gradients the energy of a mix- 
ture with doping x m i n < x < x rnax interpolates linearly 
between E m i n and E max . The effect of Coulomb interac- 
tions and gradients is to increase the total energy of these 
intermediate states to Etot{x) as shown in Fig.^. The de- 
pendence of p on x is then given simply by \i = SEtot/Sx, 
which is consistent with the numberical results of Fig.^. 

Because E to t{xmin,max) = E min , max in models with 
continuous charge one has 



(n(x) — fj, c ) dx = 0. 



(7) 



which applies beyond the SO (5) theory. Our numerical 
results are consistent with (m). 
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FIG. 4. State energy as a function of x: In the phase sep- 
aration regime, the energy becomes a convex function of x. 

Experimentally, investigations of the chemical shifts in 
LSCO |2Ul and BSCCO fM have shown that while the 
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shift is large in overdoped samples, it is strongly sup- 
pressed and pinned in underdoped samples, in agreement 
with the phase separation picture. However, due to poor 
experimental resolution, it is not possible to ascertain 
the non-monotonic behavior of fi(x). More experimen- 
tal work is needed to test this prediction of our model. 
While this work was motivated by experiment, it should 
be emphasized that extensions of our model would be 
needed to make real contact with experiments. The lat- 
tice anisotropy of the cuprates will lead to an anisotropic 
AF/SC interface energy. We expect this anisotropy to 
enlarge the region of striped phase stability relative to 
that of the droplet phases, as the former can take best 
advantage of that anisotropy. Also, disorder will make 
the coefficients J and K (or Xab and g) and the charge of 
the counterions position-dependent. Although the result- 
ing effects are complex in character, we may speculate 
that for small disorder, the defects act as pinning cen- 
ters for the stripes and lead to distortions of the domain 
structure as well as a loss of long-range order. This may 
explain the failure to observe droplet phases in the high- 
T c superconductors. For strong randomness, the size of 
the domains would be predominantly set by the disorder 
instead of the long-range interactions ^7j,^8|. However, 
we expect that the linear relation between the superfluid 
density and the doping should still hold in these glassy 
materials. 

In summary, we have shown that the interplay between 
the long-range Coulomb interaction and the local ten- 
dency to phase separation of the SO (5) model leads to 
an interesting and remarkably rich phase diagram for the 
clean system. We found that the frustrated phase separa- 
tion between hole-rich and hole-poor regions can provide 
an explanation for the gross features of the cuprates near 
the AF/SC transition when lattice anisotropy and impu- 
rity effects are taken into account. However, the SO(5) 
energy functional cannot have topological solutions as its 
lowest energy state. Therefore we believe that the topo- 
logical nature of stripes that are observed in experiment 
must arise from microscopic properties of the coexisting 
states. Finally, we draw attention to the behavior of the 
chemical potential in the phase separation regime. 
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